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An efficient method is deveioped which aiiows the calculation of distributions of the surface 
density in the equilibrium disk configurations with an isolated point mass at the center corresponding 
to known distributions of the angular velocity. We demonstrate the existence of the upper limit for 
the ratio 'mass of the galactic disk / mass of the central black hole' which essentially depends on 
the form of rotation-curves. 
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The super-massive black holes at the nuclei of galaxies and quasars are able to provide a high luminosity of 
these objects thanks to the disk accretion (Lynden-Bell 1969). Nowadays there exist serious arguments in favor of 
that hypothesis for a number of disk-like galaxies (Rees 1998; for the discussion of the observational data providing 
evidence about the existence of compact nuclei in spirals see Rubin and Graham 1987; Sofue 1996; Ratnam and Salucci 
' 2000), so the problem of finding the mass-density distributions and total masses of flattened galaxies possessing a 
central black hole from the known distributions of the matter rotation-velocity is of great astrophysical interest. The 
rotation-curves are normally constructed from the measurements of the Doppler effect for the lines 21cm, H a etc. 
(see, e.g., Schmidt 1957; de Vaucoulears 1959; Burbidge et al. 1960; Rubin et al. 1980; Carignan and Freeman 1985; 
1 Sofue and Rubin 2001; Persic and Salucci 1988). The mathematical problem of reconstruction of the surface mass- 
density from the rotation-curves in infinite disks without a central body was discussed in the papers by Burbidge 
et al. 1960; Brandt 1960; Brandt and Belton 1962; Toomre 1963; Binney and Tremaine 1987 in relation with the 
problem of non-radiating mass in the disk-like galaxies. Mestel 1963 proposed models of flattened galaxies as disks 
of finite radius with constant angular or linear rotation velocities. 

Another important physical aspect of the problem (we do not touch it in this paper) is related to the formation of 
the exterior massive parts of accretion disks in which the attraction to the disk because of its mass becomes greater 
or of the order of the vertical component of the attraction force from the part of the black hole. In this case the 
exterior part of the accretion disk considerably swells, and the proper field of the disk may play a substantial role in 
the thermal balance, distribution of pressure, etc. (Paczyhski 1978; Kozlowski et al. 1979, Kolykhalov and Sunyaev 
1980). For millions years of accretion the mass of the accretion disk can reach 0.001 the mass of the black hole due 
. to transition of the angular momentum into the exterior part of the accretion disk (Kolykhalov and Sunyaev 1980). 
In our paper the mass of the accretion disk around a black hole is neglected, and we only take into account the 
proper gravitational field of the galactic disk. In the latter we neglect the influence of viscosity, pressure and all 
kind of non-circular movements of matter, as it is done in the classical articles (Schmidt 1957; Burbidge et al. 1960; 
Brandt 1960) because of a small relative dispersion of velocities. It should be emphasized that the gravitational field 
of dark matter and that of the non-disk component (bulge) exert the same influence on the rotation-curves as the 
galactic radiating disk itself. The phenomenological expressions for these potentials are given in the book of Binney 
and Tremaine 1987 and in the paper by Lovelace et al. 1999, but some researchers are rather sceptic about the 
methods of separation of those fields (Burstein and Rubin 1985; Persic and Salucci 1988). 

The aim of our paper is to solve a new Newtonian notential-theory problem for the gravitational field of a disk 
extending at some distance from the central point mass.u The rotation-curve in the disk is supposed to be known, and 
it is necessary to find the corresponding mass of the central object and the distribution of the surface mass-density in 
the disk. We give the general solution of this problem in terms of two successive quadratures (20), (21) and formulas 
(19), (23) for the masses of the black hole and disk, respectively. As an illustration of the application of general 
formulas we derive a large class of the asymptotically Keplerian angular-velocity distributions and demonstrate a 



*On leave of absence from the Department of Hydrodynamics, Lomonosov Moscow State University, Moscow 119899, Russia 
1 Mention that the disks with an isolated central mass cannot be obtained as limiting cases of confocal spheroidal shells of 
small eccentricity, unlike the continuous disks previously considered. 
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strong dependence of the upper limit of the masses of disks on the choice of a class of rotation-curves. In spite of a 
number of strong assumptions (neglect of a non-disk component, extrapolation of data about the circular motion of 
matter into the outer parts of galactic disks where there is no radiating gas) our formulation can be considered as a 
possible approach to the solution of the problem of a hidden mass in the galactic disks with a black hole at the center. 



II. THE METHOD 

Recall that in the axisymmetric case a thin galactic infinite disk with an attracting central body of mass M can be 
described by the following Newtonian potential (in cylindrical coordinates p, z) 

°{Po)Po dp d^ GM 



*{P,*) = 2G I , 2 +— r = V^T^ (1) 
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where G is Newton's gravitational constant. The matter of the disk is distributed in the plane z = in the region 
exterior to the circumference of radius a. From the condition of balance of the gravitational and centrifugal forces 
one gets 
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(here we do not consider other forces like pressure of the gas and the radiation pressure). 

Substituting the potential (1) into Eq. (2), one easily arrives at the integral equation with the kernel divergent at 
P = Po- 

OO 

/GM 
cr (Po)K(p, p )dp = lu 2 p — (3) 
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where 

P v , s signjp- Pq )E(t) K(t) 4p p 

— K(p,p ) = +| -, T=— rj, (4) 

2po P + PQ I/" -Pol {P-Pa) 

E{t) and K (r) being the complete elliptic integrals. 

Eq. (3) is a complicated singular integral equation of the non-Fredholm type. Even the problem of finding ui for 
a given surface density a would meet difficulties during its numerical resolution. However, from the point of view of 
the observational astronomy, it is likely to solve the inverse problem, i.e., to determine the mass distribution a for 
a known to since from observations it is difficult to establish the matter distribution in the galactic disk because of 
the problem of a hidden (non-radiating) massp whereas the rotation-curves are able to provide information about 
Keplerian frequencies (Schmidt 1957; Burbidge et al. 1960; Brandt 1960; Binney and Tremaine 1987) thanks to the 
combined attracting field. Viewed in this way (finding a from uj), the straightforward resolution of Eq. (3) looks 
hopeless even by numerical means. 

To circumvent the need of solving the cumbersome Eq. (3), instead of (1) we shall take (j> in the form 



If r GM 
(ft(fi, z) = ^ / d0 / ln[(s - p cos Of + z 2 } a(s) ds + — , (5) 



where the real function a(s) has the meaning of the density of sources distribution. 

This potential satisfies Laplace's equation everywhere except in the plane z = 0, the sources of the disk being 
distributed outside the circle p < a. The main advantage of such a representation over (1) is that it allows to 



2 In practical approaches the surface density is normally reconstructed with the aid of the photometric data (Freeman 1970; 
for a critical discussion of this approach see Burstein and Rubin 1985; Persic et al. 1996). 
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introduce into the problem under consideration the powerful theory of analytic functions in the complex plane since 
the first term in (5) can be cast into the form f£ f(z + ipcos6>) d9. 
Let us show that 



Indeed, from (5) we have 

d<t> _ 1 
dz ~ 2^ 
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Tending z — > +0 for pcosO = sq > a and using the Sokhotsky-Plemelj formula for the Cauchy-type integrals, 
namelytj 



f a(s)ds f a(s)ds . 
hm / =/ — ±ma{s ), 
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we get (in our case Z — sq ±iz) 



dz 



0, 



< s < a 



\ I a(s )d0, s > a. 
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From (9) the formula (6) follows immediately. 

Let us consider the dependence of on p when p > a and z — ► 0. From (8) we obtain 
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Using now the formula 
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Eq. (10) can be rewritten in the form 
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Let us show that for a finite mass of the disk, the integral J°° a(s)ds should vanish. 
Indeed, from (6) we have for p > a 



-4nGa{p) 
hence the mass of the disk is equal to 
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3 Here so is a point on a smooth contour C: so £ C; Z — *• +so means that the point Z is approaching the point so from the 
left side in the sense of a positive running of the curve C; the symbol-^ denotes the principal value of the respective integral. 
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In order the mass of the disk be finite, it is necessary that 
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Substituting now (12) into the equilibrium condition (2), we obtain 



, GM 1 f a(s)sds 
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After introducing the new variables t — a 2 /s 2 , x = a 2 /p 2 , the last equation reduces to the Abel integral equation 

a(t) dt _ a 3 u 2 
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and the well-known formula for the solution of this equation yields (see, e.g., Sneddon 1956, p. 318) 



dt ( a 3 u> 2 



GM . 



a 2 a(x) 

2x 3 / 2 = dx J ^x~t V i 3/ 
o 

The substitution of (18) into the condition (15) yields the formula for the mass of the central body 
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hence the function a(x) can be expressed through only the known distribution of the angular velocity: 

l 
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to 2 dt 
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The corresponding distribution of the surface density, obtainable from (13), assumes the form 



(20) 
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From (21), (20) follows that the asymptotics of the surface density near the inner rim of the disk p = a is given by 
the formula 



-Att 2 Gct{x) w Ay/1 - x , A = a 



UJ 

W 2 



1 I u 2 G(M + M d ) 

2 1*372 ^3 



dt 



y/T^t 



(22) 



It can be seen from (22) that in the presence of a central body possessing a finite mass the surface density near the 
inner rim of the disk has no finite limit when a tends to zero. The inequality A < imposes a physical restriction on 
the angular-velocity distribution. 

Therefore, we have finally succeeded in expressing the mass-density distribution in the disk exclusively through the 
supposedly known distribution of the angular velocity via two quadratures which can be easily taken numerically in 
the general case, and analytically for large classes of particular distributions. 
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The total mass of the disk, Md = 2ir f up dp, is obtainable from (18) and (21): 



i * a 2 a(t) a 3 uj 2 „, , 

GM d = lim — = lim — -r- - GM . 

t _ 2t t->o t 3 / 2 

To give an analytical illustration of the results obtained, let us consider the following smooth distribution of 
angular velocity possessing the Keplerian asymptotics:El 

2 3 CL 2 Q 3 d n ^-^\ ,h/2 

uj p = a Q + a 2 ^r + a 3 ^r + . . . + a n — = a + } a k t 1 , 
P P 3 P ^ 

where a k are some constant coefficients. 

Substituting this expression into (18) we arrive at 



a(x) 



JL {ao _ GM ) + £(* + 1)x^-WbQ-, 1 + ~ 



k=2 



B(x,y) denoting the Euler beta-function. 

From the condition (19) we obtain the mass of the black hole: 



2GM = 2a Q + J2 B (l^ + ^)^ 



k=2 



Substituting the expression (25) into formula (21) one finds 
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t k / 2 dt 
y/t — X 



-47T 2 G^l=4(a -GM)VT^+J2(k + l)B(^,l + ^y k J k , J k = f 

Using integration by parts, for the integrals J k it is easy to obtain the recurrent formula 

(k + 1) J fc = 2Vl-x + kxJ k -2 , 
the use of which in (27) gives the explicit expression for the surface density: 

, , i 3/ Vl - x /l k\ ( k k(k-2) rt/21-i 

^ = ~—2^G- ]2 B (r 1+ 2h{— +X (k - l)(fc - 3) + ■ ■ • + X 



k(k - 2) . . . (k - 2[k/2] + 2) f i (k - 2[k/2])x ^ 1 + VT^x 



(fc-l)(fc-3)(fc-2[fc/2] + l) V VT^~x~ V% 

The inequality A < (cf. (22)) imposes the following restriction on the parameters a k : 



fc=2 



k \ k\ 1 \ 
-B -, - I B -, 1 

2 12 2/ 2 \2' 



< 0. 



Mention that when the function oj 2 p 3 is representable by a polynomial in inverse powers of p 2 , a2 k +i = 0, 
1,2, .:,N, then the surface density is given by the formula 



o-(ac) = - — i\ X V (3 s x s : , 

s — l 

ft = ff 2 fc (fc!) 2 (-ir +1 a 2fe 



- m=o ( 2fe ~ - s ) ! ~ m)!(2fe + 2m - 2s + 1) 



4 In (24) the coefficient ai is equal to zero by virtue of equations (17) and (18). 



III. SOME EXAMPLES 



(A) Let us consider the case when the rotation-curves belong to the family of curves with Keplerian asymptotics 
of the type 

lu 2 P 3 = GM + GM d (l + b 1 x + b 2 x 2 + b 3 x 3 ), a 2i = GM d bi, x = a 2 / p 2 . (32) 
From the formula (26) follows that 

3 4 24 

^—i-i'-s'*- < 33 > 

Taking into account (33), one obtains from (31) the distribution of the surface density corresponding to (30): 

16, 64 , /64, 128, \ 256 



. aM d Vl - a 2 p 2 

<r(p) = y — 3 

P 3 



-bi H b 3 — ( — bo H b 3 I x b-tx 2 

45 2 175 3 V 45 2 175 3 175 3 



(34) 



Formulas (32)-(34) fully describe the two-parameter family of the angular-velocity distribution and the corre- 
sponding distribution of the surface density. In Fig. 1 we have plotted the region D of the plane (b 2 ,b 3 ) on which 
the right-hand side of (34) is a positive quantity, so that (32) has physical sense. The boundary of the region D 
consists (i) of two line segments AB and AC tangent to the ellipse at the points B and C, their defining equations 
are 1 + j^b 2 + = and 1 — f^b 2 — §f&3 = 01 the P ar t °f the ellipse between the points B and C; the 

equation of the ellipse results by setting to zero the discriminant of the quadratic polynomial at the right-hand side 
of (34). 

Let us call p(b 2l bs) the minimal value of the cubic polynomial f(x): 

/3 4 24 \ 

f(x) = l-(- + -b 2 + — b 3 j x + b 2 x 2 + b 3 x 3 (35) 

on the interval (0, 1]. Then the maximally possible mass of the disk can be calculated from the condition of the 
non-negativeness of the right-hand side of (32): 

Mdmax(b 2 , 63) - -M/fi(b 2 , 63) ■ (36) 

The polynomial f(x) can assume its minimal value either inside of the interval (0, 1) or at the boundary x = 1. 
In Fig. 1 the dashed line separates the region D into two parts D\ and D 2 . Inside of Dx, where /(l) is the minimal 
value of f(x) on the interval (0, 1), the function p{b 2 ,bz) has the simple form p = — \ + jrb 2 + ^b 3 . 

In Fig. 2 the curves M drnax are shown as functions of b 2 for different fixed values of b 3 , the points {b 2l b 3 ) belonging 
to D. A simple investigation shows that the maximally possible mass of the disk can reach 9M. It is interesting 
to point out that the more complicated structure has the angular-velocity distribution in the disk (i.e. the greater 
is n in the polynomial (31)), the larger mass the disk can have in principle: for the one-parameter family of the 
angular-velocity distribution the maximal mass of the disk is equal to 5M, and for the two-parameter family the 
maximal value elevates up to 9M. 

Let us consider three particular cases of the distribution (32). _ 

1) The simplest particular case, b 2 = 63 = 0, was found by Lemos and Letelier 19940. In this case 



W = GM + GM d (l-±*\ *(p) = MMd X P ' 4 T • (37) 

V 2 p A I 7T^p 4 

From (35) and (36) follows that the mass of the disk cannot exceed two masses of the central body. 
2) The particular case b 2 — 15/8, b 3 — corresponds to the distribution of the angular velocity and surface density 
of the form 



15 a 4 \ , . 8aM d {p 2 - a 2 ) 3 / 2 

^7r a[p) = 3W 



j 2 p 3 = GM + GM d (l — 3— + -77-4 ) j °{P) = """"T 5 (T ' > ( 38 ) 



5 The inequality A < is equivalent to the condition 1 — ^62 — ||&3 > 0. 

6 See the discussion of the corresponding gravitational potential in the paper by Semerak and Zacek 2000. 
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and the mass of the disk does not exceed 5M. 
3) The choice b 2 = 45/8, b 3 = -35/16 gives 

and the mass of the disk cannot be greater than 6.69M . 

(B) To demonstrate that the mass of the galactic disk can in principle be arbitrarily large compared to the mass of 
the central black hole, let us consider the following rotation-curvesEI 

Lu 2 {p 2 +c 2 ) 3/2 = G(M + M d ) , c = const. (40) 

In this case the equation (23) is fulfilled automatically. Let us denote the parameter c 2 /a 2 as b. 
Calculating a(x) according to formula (20), we get 

a(x) _ b(-l + (3 + b) + bx 2 ) G (M + M d ) 

„3 ' \ ' 



2ax 3 / 2 ^(1 + b)(bx + l) 2 a 3 

Substituting the expression (41) into (21), we obtain the formula for the distribution of the surface density: 



(M + M d )Vbx3 ( y/b(T^x) Arctan v /b(l - x)/(l + bx) \ 
a[X) ir 2 a 2 y(l + b)(l + bx) [l + bxf/ 2 J' 1 ' 

Then from (19) follows a surprisingly simple formula for the mass of the disk: 

M d = bM . (43) 

The parameter b in formulas (41)-(43) can be an arbitrary positive number. Therefore, the mass of the galactic 
disk can be by far greater than the mass of the black hole (if the linear velocity V = puj achieves its maximum at the 
distances much larger than the accretion part of the disk 6 > 1). Mention that when a — > 0, the variable x tends 
to zero too, while the parameter b tends to infinity. From (43) then follows that M — > 0, and from formula (42) we 
obtain Kuzmin-Toomre's result for the surface-density distribution (40): 

/ , cM d 

g(p) - 2 7 r(p 2 + c 2 )3/ 2 - ™ 



IV. CONCLUSION 

Therefore, we have succeeded in solving the problem of the determination of the surface density in massive galactic 
disks with a central massive black hole from the known rotation-curves. The problem is reduced to finding two 
successive quadratures (20) and (21) which can be easily calculated by numerical means in the most general case of 
smooth distributions of the angular velocity. Besides, we have shown how the knowledge of a rotation-curve permits 
to calculate the mass of the central body and that of the galactic disk. 

The analysis of the formulas obtained gives rise to the following two general observations. 

(i) The mass of the disk around a central body in each case should be less than some upper limiting value which 
depends on the degree of concentration of matter at the inner boundary of the disk. The mass of the disk strongly 
depends on the form of the rotation-curve and in principle can exceed the mass of the black holes many times. 

(it) For an arbitrary non-degenerated smooth distribution of the angular velocity in the disk, the surface density 
at the inner boundary tends to zero like yl — a 2 p~ 2 . It corresponds to the parabolic profile of thickness of the disk 
for a finite value of the volume density at the inner boundary of the disk. At infinity the surface density decreases as 
p~\ 

Lastly, it should be remarked that the method developed in this paper is not restricted to only the case of infinite 
massive disks with attractive central body, but also permits to treat the cases of finite disks and systems of rings. 
These last cases are in the stage of analysis and preparation now, and will become a subject of future publications. 



7 The rotation-curves (40) were first considered by Kuzmin 1956; they enter as a particular case into the class of curves 
possessing the Keplerian asymptotics at p — > oo proposed by Brandt 1960 for infinite disks without a black hole at the center. 
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FIG. 1. The region of positiveness of a(p) in the plane (62,63). 
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